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Novel ferromagnetism in Cai-aiLa^^Be is studied in terms of the Ginzburg-Landau theory for ex- 
citonic order parameters, taking into account symmetry of the wavefunctions. We found that the 
minima of the free energy break both inversion and time-reversal symmetries, while the product 
of these two remains preserved. This explains various novelties of the ferromagnetism and pre- 
dicts a number of magnetic properties, including the magnetoelectric effect, which can be tested 
experimentally. 



Novel ferromagnetism in Cai-ajLaj^Bg ||l| has been the 
subject of extensive studies due to its high Curie temper- 
ature (~ 600K) in spite of a small moment (0.07/iB/La) 
and lack of partially filled d- or /- bands. However, its 
mechanism still remains controversial. Ceperly suggested 
that this is a first example of the ferromagnetic phase of 
a dilute electron gas . An alternative explanation [g[Q] 
is based on the excitonic state g of the parent material 
CaBg. The parent compound CaBg has a cubic struc- 
ture, and some band structure calculations ||,|j predict 
a small overlap of the conduction and valence bands at 
the X-points. Matrix elements of the dipole moments 
vanish between these states at the X-points; the dielec- 
tric constant is not enhanced even when the band gap 
collapses. This makes the excitonic insulator a plausible 
candidate for the ground state of CaBg . By La-doping 
the extra electrons are doped into this excitonic insulator, 
and it has been found in the mean-field approximation 
that these electrons are perfectly spin-polarized |^ . Nev- 
ertheless, there is so far no conclusive explanation for the 
magnitude of the moment, much smaller than the elec- 
tronic moment doped by La. Furthermore, the correla- 
tion between the La-doping concentration and the ferro- 
magnetic moments are questioned by recent experiments 
g . Several experiments suggest that the ferromagnetism 
is not a bulk phenomenon, and occurs only in the thin 
film sample |9|, or near the surface as evidenced by an 
electron spin resonance (ESR) experiment |lO| . 

In this letter we study in depth the symmetry prop- 
erties of the excitonic state in CaBg in terms of the 
Ginzburg-Landau (GL) theory, and propose a possible 
scenario for the novel ferromagnetism. We classify pos- 
sible states in terms of the magnetic point groups. The 
idea is that CaBg is a triplet excitonic insulator [|| with 
broken time-reversal (i?) and inversion (/) symmetries, 
while their product RI is kept intact g]. This means 
that CaBg is an antiferromagnet (AF), and the ferro- 
magnetism is induced by the magneto-electric (ME) ef- 
fect, and the La-doping and/or the surface works mainly 
as a source of the electric field. We emphasize that the 
present theory of ferromagnetism is different from those 



in l|,Q], though all these theories are based on exciton 
condensation. 

This scenario of exciton condensation is justified when 
the band structure has a small overlap/gap at the X 
points. It is, however, still in controversy whether it is 
the case. Contrary to the LDA calculation |§,0|, more re- 
cent LDA-I-GW calculation [|l| suggested a large gap of 
0.8eV, too large compared with the energy scale for exci- 
ton binding energy (~ 600K). An angle-resolved photoe- 
mission spectrum (ARPES) ||l2[ also shows a large gap 
of order leV, and is in good accordance with the GW 
result. Nevertheless, there are still enough evidences for 
believing a small gap or a small overlap. First, another 
recent GW calculation [Q shows a small overlap at the 
X points. Second, de Haas-van Alphen measurement |1J] 
finds Fermi surfaces of both electrons and holes. Third, 
the ARPES experiment is surface-sensitive, and one can 
argue that the ARPES result |l2| of the large gap might 
be attributed to surface effect. Fourth, measurements 
of X-ray scattering and Raman scattering in CaBg Il3] 
shows an anomaly at 600K, below which the compound 
is determined to be tetragonal. This indicates that the 
material has some ordering and symmetry-breaking even 
in the parent compound. This ordering is antiferromag- 
netic, as indicated in /iSR measurements [|l6| . Therefore, 
a large gap is unlikely, which prohibits any instability 
and sensitivity to the small amount of doping. Thus, al- 
though there is still no concensus on a size of the gap, the 
excitonic insulator still remains one of the most promis- 
ing candidates for the novel ferromagnetism. 

Let Qa;, Qy, Q^ denote the three X points 
(7r,0, 0), (0,7r,0), (0,0,7r) in the cubic Brillouin zone. 
When excitons are formed, excitonic order parameters 
(^ka'^k'/?) will have nonzero values, where 6kcr and a^^ 
are annihilation operators of electrons with spin a at the 
conduction and the valence bands respectively Q. Be- 
cause the excitonic instability occurs only in the vicinity 
of the three X points, we make the following assumptions; 

(i) The order parameters are k-independent near the 
X points. 



(ii) The order parameters connecting different X points 
are neglected Q . 

As a result we keep only the order parameters //(Qi)^^ = 

i^Qic^QiP)- The cubic (Oh) symmetry of CaBg restricts 
a form of the GL free energy $, as is similar to the GL 
theory for unconventional superconductors |jl^,|l8| . 

At the X points, the k-group has a tetragonal (D4h) 
symmetry, and the conduction and the valence band 
states belong to X3 and Xg representations, respectively 
1^,0. Without the spin-orbit coupling, the order param- 
eter ri{Qi)ai3 belongs to X3 x Xg = X'^. From now on we 
shall take the spin-orbit coupling into account. Then, the 
point-group transformation is accompanied by the spin 
rotation, and the representation of rj will be altered. 

Let us consider the triplet channel for the excitons, be- 
cause the exchange interaction usually favors the triplet 
compared with the singlet [^. Then the spin-1 represen- 
tation of the rotation group will be multiplied, and we get 
X4 + X5 in the D4h group. Explicitly, for the triplet ex- 
citons at the Q^, the Sz = component rj^i — f^n obeys 
the X4 representation, while the Sz — H components 
rjil, ?7|| obey the Xg. Here the spin-quantization axis is 
taken to be -l-^-direction. Since we are taking into ac- 
count the spin-orbit coupling, these up- and down-spins 
should be interpreted as pseudospins p^ ]. 

When we consider the three X points, the cubic sym- 
metry is restored since the order parameters have a 
wavevector q = 0. The order parameters have nine 
components in total, which can be classified into irre- 
ducible representations of Oh- They split into three 3- 
dimensional representations: Fis + Fis + F25. Let us 
call basis functions as 77i(Fi5,l), r7i(Fi5,2) and r7j(F25), 
where i = x,y, z and j = x{y'^ — z'^) , y{z'^ — x'^) , z{x^ — y^). 
These basis functions transform like a vector (a;, y, z) for 
Fi5 and like {x{y^ -z'^),y{z'^-x^),z{x'^-y'^)) forFzg. By 
judicious choice of phase for basis functions, they can be 
made to transform under the time-reversal R as complex 
conjugation, and they are odd under the inversion /; i.e. 
It] = — ry, and Rq = r\* . The basis functions are given as 



?7z(ri5, 1) = — =(77ii(Q2) - ??tt(Q^))> 

r;,(Fi5,2) = -i{i(?7Ti(Qy) + 7?iT(QJ) 
?7^(x2-y2)(F25) = 2{i('7n(Qa)+^iT(Qa)) 



(1) 



(2) 



(3) 



Here, the spin-quantization axis in ?7Q,/3(Qi) is taken as 
-fi-axis (i = x, j/, z). Other components are obtained by 
cyclic permutation of x, y, z. 

Let us now write down the GL free energy in terms of 
these order parameters. The GL free energy $ should be 
invariant under the elements of Oh and under the time- 
reversal R. We make two remarks helpful in writing down 



$. First, only even-order terms in r/ are allowed by the 
inversion symmetry. Second, owing to the time-reversal 
symmetry, the order of Imr/ in each term should be even. 
Thus, $ is given up to quadratic order as 

$(2) = Ai ^(Re ,7,(Fi5, 1))' + A2 ^(Re m{V^^,2)f 

i i 

-I- ^3 ^ Re ?7»(ri5, l)Re ?y,(Fi5, 2) -I- A^ ^(Re r/, (F25))' 
+ Bi ;^(Im 7?,(Fi5, \)f + S2 ^(Im r?.(Fi5, 2))2 

i % 

+ B3;^Im 7y,(Fi5,l)Im 77,(Fi5, 2) -f S4 ^(Im 77,(F25))' (4) 

Therefore, as the temperature is lowered, one of the fol- 
lowing states will be realized; 

(A) Re r7(Fi5, 1) = cRe J7(ri5, 2) ^ 0, 

(B) Re r7(F25) ^ 0, 

(C) Im r/(Fi5, 1) = dm r/(Fi5, 2) ^ 0, 

(D) Im r7(F25) ^ 0, 

where c is a nonzero constant. Therefore, a condensation 
of excitons in F15 and in F25 do not occur simultaneously. 
In each case among (A)-(D), all directions of the vector 
7] are degenerate, and this degeneracy will be lifted in 
the quartic order, as we shall see later. All these states 
are accompanied by a lattice distortion. This distortion, 
however, is expected be rather small because it couples 
to the order parameter in the quadratic order, not linear 
order. Indeed the tetragonal distortion detected by X-ray 
scattering is as small as 0.03% [ p^ . 

The states (A),(B) preserve time- reversal symmetry 
and thus nonmagnetic. From symmetry consideration, 
this implies that neither the ME nor the piezomagnetic 
(PM) effect will be observed [Q . Roughly speaking these 
states are far from showing ferromagnetism. 

To study which state is realized, we should know the 
coefficients Ai^ Bi in eq.(Q), which are related to the fol- 
lowing matrix elements for the exchange interaction by 
the Hartree-Fock approximation ||] . 

I dx I dxV^Q. (x)04. (x) 7^^<^q1 i^W^. (x') , (5) 

where 0q . (x) is Bloch wave function of the c-band (c = 
a,b) at k = Qi, e is a dielectric constant. In the cal- 
culation of these matrix elements, we neglect the spin- 
orbit coupling. To calculate 0q. (x) we proceed as fol- 
lows. We define 4>q^X' ^^'^ 'f'o^x' ^^ linear combina- 
tions of p-orbitals of B with Xg and Xg, respectively, 
and c!)R^'y' as a linear combination of d-orbitals of Ca 
with Xg. Then, (/)q. (x) is defined as a bonding orbital of 
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q'^-j^, and ^q'J X' ' ^^^ '/'cji (x) ^^ equal to (pafxa 1 '^i^h 



proper normalization ||]. We evaluated the coefficients 
in (0), and found that even if coupling terms between 
the three X-points are included, the coefficients for the 
imaginary parts of the order parameters are smaller than 
those for the real parts . Therefore, the states (C) and 
(D) would be favorable than (A) or (B). Thus, from now 
on we shall concentrate on (C) and (D) pOJ ]. 

The states (C) and (D) break both the inversion and 
the time-reversal symmetries; they are magnetic states. 
A crucial observation is that they preserve the RI sym- 
metry PI , leading to many interesting consequences listed 
below. (I) The RI symmetry prohibits a uniform mag- 
netic moment. Dzyaloshinskii used this symmetry to ex- 
plain why weak ferromagnetism is present in a-Fe2 03 
while not in CraOa |2^. Thus the states (C) and (D) 
are antiferromagnetic. This agrees with the result of the 
/xSR measurement |16[| with a moment of 0.0039/iB/uiol. 
Note that the magnetic unit cell is identical with the 
original unit cell. Thus, no extra Bragg spots appear be- 
low the AF transition. (II) The RI symmetry also pro- 
hibits the PM effect. This follows because a stress cannot 
break the RI symmetry which impedes ferromagnetism. 
Nevertheless, a gradient of stress can break this symme- 
try and will induce ferromagnetism. We note in passing 
that if intervalley excitons condense, i.e. the assumption 
(ii) is violated, the PM effect can occur. (Ill) The RI- 
invariance results in the linear ME effect, as observed in 
Cr203 |22]. This occurs because an external electric field 
E breaks this RI symmetry and enables ferromagnetism 
p2[ . An electric field E belongs to the Pis, and couples 
linearly with the order parameters in Pis as 



5<P 



-CiE • Re r7(Pi5, 1) - C2E • Re r,(Pi5, 2). (6) 



Here the imaginary parts of the order parameters are ab- 
sent due to invariance of (5$ under time-reversal. Thus, 
in the presence of E, both the real and imaginary parts of 
the order parameters acquire nonvanishing values, result- 
ing in a uniform moment. As for the P25, similar effect 
can be found. (IV) The optical non-reciprocal (NR) effect 
in reflection is predicted to occur, as is similar to Cr203 
|P3| . This effect will occur only below the Neel temper- 
ature. Generally, i?/-invariant materials with broken R- 
and /-symmetry will exhibit the NR effect |g^ . 

Let us consider quartic order terms $'^^'' in the GL free 
energy, which lifts the degeneracy in the direction of rj. 
We do not write down its lengthy formula here |2J]. By 
minimizing $'^' -|- <&'■'*% we find four possibilities; 

(C-1) Im J7(Pi5, 1) - cilm r7(Pi5, 2) = (0, 0, C2), 

(C-2) Im J7(Pi5, 1) = cilm r7(Pi5, 2) = (c2, C2, C2), 

(D-1) ImJ7(P25)-(0,0,c), 

(D-2) Imr7(P25) = (c,c,c), 

where c's are constants. The direction of the lattice dis- 
tortion is tetragonal in (C-1) and (D-1) and is trigonal 



in (C-2) and (D-2). Magnetic point groups G for these 
states are (C-1) 4/m'mm, (C-2) 3'm, (D-1) 47m'm'm, 
(D-2) 3'm'. In each case among (C-l)-(D-2), there are 
two types of degenerate AF domains, when we fix the 
axis of tetragonal or trigonal distortion. Since there are 
three and four choices of axes for tetragonal and trigonal 
cases, respectively, total degeneracy is six in (C-1)(D- 
1) and eight in (C-2)(D-2), which is equal to an order 
of the quotient group (Oh x {E, R})/G. We can draw 
some analogies with anisotropic superconductivity (SC). 
The order parameters t] of triplet excitons correspond to 
the d-vector in triplet SC. It is nevertheless misleading 
to look for SC counterparts of our phases (C-l)-(D-2), 
because our order parameters are confined in the neigh- 
borhood of the three X points. They are triplet and even 
functions in k, which never occurs in the SC. 

Information for an ME effect can be obtained from 
ref. [|l9[. When we write a bilinear term of H and E in 
the free energy as aijHiEj, the property tensor a.y is 



ai 

(C-1) ( -ai 





(D-1) 




ai — Qfi 

(C-2) ( -ai ai 

Qfl —Qfi 

ai a2 a2 

(D-2) a2 OLx a2 

0.2 oi2 a\ 



These explicit forms of ol contain information on a di- 
rection of magnetization M under an electric field E and 
that of polarization P under a magnetic field B. In par- 
ticular, in (C-1) and (C-2), M_LE and P_LB always hold. 
Thus, measurement of the ME effect with a single crys- 
tal of CaBg will manifest which of these four is realized. 
Experimentally, domain structure of the AF develops, 
and the sign of a. is reversed when the staggered mag- 
netization is reversed. In measurement of the ME effect, 
the domain structure can be aligned by a magnetoelectric 
annealing, in which the sample is cooled under both elec- 
tric and magnetic fields. Domain boundaries between the 
two AF domains can exhibit interesting properties. As in 
boundaries between two SC domains with broken time- 
reversal symmetry |18| , localized current and magnetic 
moment are induced near the boundary. In the present 
case it is interpreted as the ME effect. 

We remark on the recent results of X-ray scattering 
and Raman scattering p5[ , which strongly support our 
theory. They show a tetragonal distortion below 600K, 
which indicates that (C-1) or (D-1) is the case. Further- 
more, inversion-symmetry breaking does not appear in 
the Raman spectrum p5| , which is consistent with (C) 
or (D), but not with (A) or (B). Thus, (C-1) and (D-1) 
are the only possibilities totally consistent with |l5[ . 

The ferromagnetism in the thin-film CaBg is inter- 
preted as caused by an electric field between vacuum and 
the substrate. For the powder experiment and the La- 
doping experiment, the explanation is more delicate. We 



believe that the carriers by La-doping are trapped by 
impurities/defects and create local electric fields. There- 
fore, an internal electric field and/or a gradient of a strain 
has a random direction, and hence the magnetic mo- 
ment is induced locally due to this mechanism. At first 
sight they appear to cancel with each other, giving zero 
or quite small uniform magnetization, which contradicts 
with the experiments. However the random direction of 
the frozen electric field and the external magnetic field 
determine the domain structure of the AF; the free en- 
ergy to be minimized is composed of (i) the energy gain 
due to the magnetization in the external magnetic field 
and (ii) the elastic energy loss of the spatial change of the 
order parameter. Therefore, the pinning of the domain- 
wall motion leads to the hysteresis behavior, which is 
regarded as the experimental signature of the "ferromag- 
netism" . The detailed study of the magnetization curve 
including its dependence on the sweeping time, however, 
can distinguish it from the real ferromagnetism. 

Other peculiarities of Cai_a;Laa;B6 can also be ex- 
plained as well. The high Curie temperature (~ 600K) is 
nothing but a Neel temperature of the parent compound 
CaBg, and is not contradictory with a tiny magnetic mo- 
ment. A rather narrow range {x < 0.01) of La-doping 
allowing ferromagnetism is attributed to fragility of ex- 
citonic order by a small amount of impurities p6| , p7[ . 
Moreover, our scenario is also consistent with the experi- 
mental results that deficiency in Ca sites psl or doping of 
divalent elements like Ba m or Sr |29] induces ferromag- 
netism. It is also confirmed numerically by a supercell 
approach that imperfections and surfaces can induce a 
local moment [ pO[ . It is hard to explain them within the 
spin-doping scenario 0. Furthermore, strangely enough, 
it is hard to find a correlation between magnetism and 
electrical resistivity, as seen in magnetizaion ^8| and in 
nuclear magnetic resonance |3^. This novelty is a nat- 
ural consequence of our scenario; electrical resistivity is 
mainly due to doped carriers, while the magnetization is 
due to local lattice distortion and/or electric field. 

The ESR experiments by Kunii [10| also support the 
above scenario. The ESR data show that in a disk- 
shaped Cao.ggsLao.oosBg, the magnetic moment only ex- 
ists within the surface layer of ~ 1.5/im thick. Further- 
more, the moment M does not orient in the direction of 
H, i.e. it feels strong magnetic anisotropy to keep the 
moment within the disk plane. This might be due to the 
long-range dipolar energy, and not due to the above sce- 
nario. Nevertheless, it is unlikely that the long-range 
dipolar energy causes such a strong anisotropy. This 
point requires further investigation. 

Let us, for the moment, assume that this strong 
anisotropy is mainly caused by the exciton condensation 
and the ME mechanism. Since this electric field should 
be perpendicular to the plane, the strong easy-plane 
anisotropy parallel to the surface implies that M_LE. 
Therefore, among the four cases, (C-1) or (C-2) are com- 



patible. Together with the result of Raman scattering, 
(C-1) is the most promising candidate for CaBg. 

We mention here a role of the spin-orbit coupling. In 
the limit of zero spin-orbit coupling, the states (C) and 
(D) become degenerate, and the quartic-order terms in 
the GL free energy do not lift this degeneracy. Sixth- 
order terms will lift it, and a resulting state will belong 
to either 4/m'mm or 4'/m'm'm. Both of them still lead 
to the ME effect in the absence of the spin-orbit coupling; 
this ME effect must be generated from an orbital motion. 
Thus, the AF state in CaBg has an orbital nature as well 
as a spin nature. With spin-orbit interaction, these two 
are inseparably mixed together. 

In conclusion, we have studied the symmetry proper- 
ties of the excitonic state in CaBg, and found that the 
triplet excitonic state with broken time-reversal and in- 
version symmetries offers a natural explanation, in terms 
of the ME effect, for the novel ferromagnetism emerging 
in La-doping or thin-film fabrication. This scenario can 
be tested experimentally by measurements of the ME ef- 
fect and the optical non-reciprocal effect in single crystal 
of the parent compound CaBg. 
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